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ABSTRACT 

We investigate the effects of the environment on star-formation in a sample of massive luminous and ultra- 
luminous infrared galaxies (LIRGs and ULIRGs) with S(24 ^m)> 80/iJy and /+ < 24 in the COSMOS field. 
We exploit the accurate photometric redshifts in COSMOS to characterize the galaxy environment and study 
the evolution of the fraction of LIRGs and ULIRGs in different environments in the redshift range z = 0.3 -^ L2 
and in bins of stellar mass. We find that the environment plays a role in the star formation processes and 
evolution of LIRGs and ULIRGs. We find an overall increase of the ULIRG+LIRG fraction in an optically- 
selected sample with increasing redshift, as expected from the evolution of the star formation rate density. We 
find that the ULIRG+LIRG fraction decreases with increasing density up to z '-^ 1, and that the dependence 
on density flattens with increasing redshift. We do not observe the reversal of the star-formation rate density 
relation up to z = 1 in massive LIRGs and ULIRGs, suggesting that such reversal might occur at higher redshift 
in this infrared luminosity range. 

Subject headings: galaxies: general -galaxies: evolution - large-scale structure of universe - surveys 



1. INTRODUCTION 

Galaxy properties such as morphology, color, star- 
formation rate do vary with redshift, stellar mass and environ- 
ment (Kauffmann et al. 2004, Postman et al. 2005, Cucciati 
et al. 2006, Cooper et al. 2008). In the local Universe, large 
redshift surveys (SDSS) have shown that the star-formation 
properties of galaxies are dependent on the local environment 
both in clusters and in low density environments, such as 
cluster outskirts and the field. Red passive early-type galax- 
ies (E+SO) are more frequently found in high density envi- 
ronments, while blue late-type, star-forming galaxies tend 
to dominate in low density regions (Dressier 1980). The 
average star-formation rate (SFR) decreases with increasing 
density (SFR-density relation), suggesting a scenario where 
galaxy interactions, and in general galaxy environment, play 
a major role in regulating star-formation. AGN activity is 
also linked to the environment. In the local Universe lumi- 
nous AGN are found preferentially in the field and in under- 
dense regions, where their hosts can preserve large amounts 
of gas to feed the black-hole (Kauffmann et al. 2004). Us- 
ing the SDSS data, Hwang et al. (2010) found no depen- 
dence of the fraction of LIRGs plus ULIRGs among infrared 
galaxies versus galaxy density at fixed stellar mass. Con- 
versely, the star formation activity of LIRGs and ULIRGs 
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are found to strongly depend on the morphology and the dis- 
tance to the nearest neighbor galaxy, indicating that galaxy- 
galaxy interactions/merging play a critical role in triggering 
the star formation activity of local LIRGs and ULIRGs. The 
morphology-density and color-density relations are observed 
at intermediate redshift and apparently are in place by z ^ 1 
(Cassata et al. 2007, Capak et al. 2007, Tasca et al. 2009, 
Postman et al. 2005, Cucciati et al. 2006). Recently, Scov- 
ille et al. (2007b) found correlations between environment 
and stellar mass, SFR, SED type and morphology for a large 
mass-selected sample over the redshift range z = 0. 1 ^ 1 in the 
COSMOS survey (Scovifle et al. 2007a). Kovac et al. (2009) 
showed that star-formation and color transformation rates of 
galaxies are higher in the group environment than outside 
groups by a factor of 2 to 4 since z ^ L Elbaz et al. (2007) and 
Cooper et al. (2008) observed the reversal of the SFR-density 
relation at z ~ 1: the average SFR increases with increas- 
ing density. An enhancement of the star-forming fraction at 
intermediate densities has also been reported for a z ^ 0.8 
cluster by Koyama et al. (2008). Tran et al. (2010) found 
an increasing fraction of IR luminous galaxies with increas- 
ing galaxy density in cluster at redshift ^L6. Caputi et al. 
(2009) analyzed the close environment, on 1 Mpc scales, of 
luminous infrared galaxies (LIRGs, L/«= 10" H- lO'^Lo) and 
ultra-luminous infrared galaxies (ULIRGs, L/r> lO'^L©) in 
zCOSMOS (Lilly et al. 2007), finding that luminous infrared 
galaxies between z = 0.6 -^ 1 are more often found in over- 
dense environments, while ULIRGs prefer underdense envi- 
ronments. Ideue et al. (2009) found an increasing fraction of 
[Oil] emitters in the high density environments at z ^ 1.2 in 
the COSMOS field, based on naiTow band imaging. 

All these results suggest that galaxy evolution and star- 
formation activity are dependent on the environment, and at 
z> 1 galaxies are caught in their active star-forming evolu- 
tionary phase in high density regions such as the center of 
groups. This phase might precede the formation of red ellipti- 
cals observed in the cluster centers in the local Universe. Both 
the SFR and the color-density relations exhibit strong stel- 
lar mass trends (Kauffmann et al. 2004, Balogh et al. 2004, 
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Haines et al. 2007). Therefore, in order to isolate environ- 
ment from mass-driven evolution (downsizing), it is crucial 
to study galaxy properties in bins of mass. 

Galaxy interactions, fly-bys and minor mergers are among 
the mechanisms expected to have an impact on evolution and 
star-formation activity of galaxies. Ram pressure stripping 
by the intra-cluster medium (ICM) might also contribute to 
regulate star-formation processes in clusters (Kapferer et al. 
2008a,b). However observations of intermediate redshift clus- 
ters tend to exclude ram pression stripping as a mechanism of 
star-formation quenching (Demarco et al.). Disentangling the 
effect of the cosmological decline of the SFR from environ- 
mental effects is not straightforward (Lilly et al. 1996, Madau 
et al. 1996, Le Floc'h et al. 2005, Hopkins et al. 2006, Caputi 
et al. 2007). Large unbiased samples are needed to isolate the 
global evolution of the SFR and mass-driven effect (Bell et 
al. 2004) from environment-driven effects. Deep, large-area 
surveys, such as COSMOS (Scoville et al. 2007a), are useful 
to address this topic, as they can sample large volumes of Uni- 
verse from moderate to high redshifts. In particular, the COS- 
MOS large area (2 square degrees in the sky, corresponding 
to physical scales of ^ 40 Mpc at z = 1) allows probing, from 
intermediate to high redshift, a wide range of stellar masses 
and galaxy densities, from the dense clusters and groups to 
the field. This is usually not achieved with deep pencil beam 
surveys. However, for such large area a complete deep spec- 
troscopic surey cannot be easily achieved. For this reason the 
COSMOS collaboration acquired deep and accurate photom- 
etry in (so far) 34 bands from far-U V to mid-infrared, includ- 
ing imaging in 16 intermediate and narrow band filters from 
the Subaru Telescope (COSMOS-21, Taniguchi et al. 2007, 
Capak et al. 2009). The survey has, on average, a photometric 
point every 200 A between the U and z bands (29 photomet- 
ric bands in this wavelength range), simulating low resolution 
spectroscopy. Thus, highly reliable photometric redshifts are 
easily computed (Ilbert et al. 2009a, Salvato et al. 2009). As 
demonstrated by Scoville et al. (2007b), at redshift below 1.2 
the quality of the photometric redshifts is sufficient to achieve 
the accuracy needed for this work, i.e. tracing the evolution of 
galaxies, AGN and star-formation processes in the large scale 
structure. 

In this paper we use the COSMOS dataset to characterizing 
the environment of bright star-forming galaxies , and study- 
ing their properties as a function of the environment. The 
plan of the paper is as follows: Sect. 2 presents the multi- 
wavelength dataset. Sect. 3 presents the sample selection, 
while Sect. 4 presents the method for density computation. 
Sect. 5 presents the galaxy properties as a function of the en- 
vironment, and Sect. 6 presents a summary of the results. A 
standard cosmology (flm = 0.3 and fix = 0.7) with Hq = 70 km 
s"' Mpc"' is used throughout. Magnitudes are given in the 
AB system. 

2. DATA 

The COSMOS field has deep multi-wavelength photomet- 
ric coverage over an area of 2 square degrees from X-rays 
to the radio. The photometric catalog includes data from 34 
broad, intermediate and narrow bands from 1550 A to 8 /im 
,and it is complete down to /"^ = 26.2 (Sa for point a point 
source, Capak et al. 2009). We use the photometric redshifts 
drawn from the catalogs of Ilbert et al. (2009a). These have 
an accuracy as good as (7^/(1 +z) = 0.012 up to /^ =24 and 
z< 1 .25. For fainter t magnitudes the accuracy of photomet- 
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Fig . 1 . — LiR versus redshift for the 24 fim sample with S(24 /im )> i 
Solid lines show the luminosity and redshift selections. 
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ric redshifts rapidly deteriorates (Ilbert et al. 2009a). This pa- 
per is limited to galaxies with /""" < 24, i.e. the sample with the 
best photo-z accuracy available. Ilbert et al. (2009a) photo- 
z were complemented by those computed by Salvato et al. 
(2009) for the AGN detected in the XMM-Newton survey of 
COSMOS (Cappelluti et al. 2009, Brusa et al. 2007). The lat- 
ter account for AGN time-variability and use a library of com- 
posite template SEDs, including both AGN and galaxy contri- 
butions. The precision of AGN photo-z is ct,/(1 +z) = 0.014 
for /"^ < 24 at any redshift. 

The stellar masses were derived from K-band luminosity 
and mass-to-light ratio using the Arnouts et al. (2007) re- 
lation and Salpeter IMF (Ilbert et al. 2009b). The measured 
mass-to-light ratios are consistent within 0. 1 dex with the an- 
alytical relation of Arnouts et al. (2007) for stellar masses 
larger than 3 x 10 '"M0 for both quiescent and star-forming 
galaxies (Ilbert et al. 2009b). 

2.1. MIPS 24 iim data 

The full 2 deg2 COSMOS field was observed by 
Spitzer/MIPS down to flux limits of 80 ^Jy at 24 fim (S- 
COSMOS, Sanders et al. 2007). The multiband catalog of 
the counterparts of the MIPS 24 ^m sources is described in 
Le Floc'h et al. (2009). Here we briefly recall its main fea- 
tures. The counterparts of the 24 ^m sources were first iden- 
tified by cross-matching with the /Tj-band COSMOS catalog 
(McCracken et al. 2010) using a matching radius of 2". The 
identification was based on the closest TTj-band source to the 
MIPS position, and using the IRAC 3.6 /im as a prior in case 
of multiple possible counterparts in the /Ts-band. About 90% 
of the MIPS sources above 80 /iJy have been assigned a se- 
cure near-IR counterpart, and about 5.5% have less secure 
identifications. The final catalog consist of 39360 objects. 
The /Ts-band sources were merged with the photometric red- 
shift catalog of Ilbert et al. (2009a) and Salvato et al. (2009) 
with a matching radius of 1", in order to associate a photo- 
metric redshift to each 24 /im source. Here again the closest 
/"^-band source was chosen as the counterpart of the /Tj-band 
source, yielding 90% of optical identification for the full 24 
/im catalog down to t = 26.5 (see Le Floc'h et al. (2009) for 



Environment of LIRGs in COSMOS 




2000 



1500 



g 1000 

3 



500 



X-ray detections 





Fig. 2. — Upper panel: fraction of optically-selected galaxies with mag- 
nitude I*' > 24 as a function of the stellar mass and redshift. Redshift is 
increasing from z = 0.4 to 1.2, from left to right. Lower panel: fraction of 24 
^mselected galaxies with S(24/im) > 80 fj,]y and i^ > 24 as a function of the 
stellar mass and redshift. 



a detailed discussion). The MIPS 24 /im sample used in this 
paper includes AGN detected by XMM-Newton (about 1200 
with L2-iokgv> 10"^^ erg/s, Cappelluti et al. 2009, Brusa et al. 
2007). 

Overall, the photometric redshift accuracy for the 24 /im 
sample is similar to that of the optical sample, <t~/(1+z) = 
0.010 for /+ < 24 (Ilbert et al. 2009a). In particular the frac- 
tion of catastrophic outliers is < 1% for sources with flux 
S(24fim) > 80 /iJy. 

2.2. Star-formation rates 

Chary & Elbaz (2001) have shown that the mid-infrared 
flux is a good estimator for the total infrared luminosity (L/r) 
emitted by a galaxy in the range 8 ^ 1000 /im in the local 
universe. Our sample is limited to z < 1 .2, for which the mid- 
infrared and radio estimators provide consistent L/r values at 
least in the LIRG regime (Elbaz et al. 2002). This assumption 
is based on infrared spectral energy distributions (SED) typi- 
cal of luminous star-forming galaxies, whose bolometric out- 
put is dominated by star-formation. SFR of sources detected 
in X-rays may be overestimated due to the AGN contribution 
to their infrared SED. 

Based on this assumption, we use the mid-infrared flux at 
24 /im and Chary & Elbaz (2001) templates to compute the 
L/R. This method uses a library of 105 templates and the 
24 /tm flux density to compute the mid-infrared rest-frame 
luminosity, and derive L/r for each galaxy. Figure [T] shows 
the L/R as a function of the redshift for the MIPS 24 /im 
catalog down to 80 /iJy. We convert Ljr into the total star- 
formation rate (SFRjn) by using the Kennicutt (1998) relation: 
SFRiR = 1.72 X 10-'°L/r[Lq]. We neglect the UV contribu- 
tion since this makes up a small fraction (about 10%) of the 
bolometric luminosity for LIRGs at z '--^ 1 (Burgarella et al. 
2007). The SFR are derived assuming a Salpeter IMF, consis- 
tent with that used to derive steUar masses. 

3. SAMPLE SELECTION 



Fig. 3. — Open histogram: all sources with t <24, M« > lO'^-^M©; filled 
histogram: LIRGs with the same selection; small panel: LIRGs detected in 
X-rays by XMM-Newton. 

The following analysis is based on a mass-selected samples 
drawn from the MIPS 24 /im and optical catalogues (Sanders 
et al. 2007, Le Floc'h et al. 2009, Ilbert et al. 2009a) in the 
z = 0.3 -^ 1.2 redshift range. In order to assess the complete- 
ness of our samples in stellar mass, we compute the fraction 
of optically and 24 /im selected galaxies with apparent mag- 
nitude r larger than a given value as a function of their stellar 
mass. Figure |2] upper panel, shows the fraction of optically 
selected galaxies with magnitude i^ > 24 as a function of the 
stellar mass and redshift, in four redshift bins from z = 0.3 
to 1.2. For M* > 3 X lO'^M© the fraction of galaxies with 
t > 24 is < 0.01,0.03,0.22 and 0.4 at z = 0.3 4-0.6,0.6 ^0.8, 
0.80-^ 1 and 1 ^ 1.2, respectively (Table 1). Therefore for 
/+ < 24 the stellar mass completeness at M* = 3 x 10 '^M© 
is 99, 97, 78 and 60% in each of the previous redshift bins, 
respectively. The 24 /tm detected galaxies are a subsample 
of the previous one, since they were selected to have an op- 
tical counterpart brighter than /"^ = 24. Figure |2l lower panel, 
shows the same fraction as before for galaxies with S(24/im 
)> 80 /iJy and magnitude /"^ > 24. The fraction of such galax- 
ies withM* > 3 x IO'^Mq is < 0.001, < 0.1,0.06 and 0.25 at 
z = 0.5,0.7,0.9 and 1.1, respectively. Thus, the completeness 
in stellar mass is larger than 90% up to z = 1 and 75% up to 
z=1.2. 

The cut in stellar mass adopted here implies that all star- 
forming galaxies in the sample are LIRGs or ULIRGs, given 
the relation between SFR and stellar mass at z '^ 1 (Elbaz et 
al. 2007). We investigate the effects of environment on dusty, 
star-forming galaxies with t < 24, L/r> IO^Lq and stel- 
lar mass M* > 3 x 1O'"M0 up to z = 1. In the redshift bin 
z = 1 ^ 1.2 we will consider galaxies with L/r> 2 x 10"Lq 
(Figure [T]). This selection implies a SFR sensitivity SFR> 
18M0/yr and 35M0/yr up to z = 1.2. The selected sources 
are LIRGs and ULIRGs. ULIRGs are a small fraction of the 
whole sample (^ 1.8%), therefore we will refer to the sam- 
ple as LIRGs throughout the paper, keeping in mind that it 
includes also a small fraction of ULIRGs. The properties 
of LIRGs will be compared to those of optically selected 
galaxies in the same range of stellar mass and magnitude 



Feruglio et al. 



TABLE 1 

Stellar mass completeness. 




Redshift 0.3^0.6 0.6^0.8 


0.8^1 


1^1.2 


Fo;„(> 1O'"M0) 0.06(94) 0.13(87) 
Fop,{> 3 X lOi°M0) 0.006 (99.4) 0.03 (97) 
Fir{> 3 X 1O'°M0) 0.001 (99.9) 0.01 (99) 


0.34 (66) 
0.22 (78) 
0.06 (94) 


0.57 (43) 
0.41 (59) 
0.25 (75) 



Note. — Fraction of galaxies in the optical and 24 fim samples with ('+ > 24 and stellar mass larger than lO'" and 3 X lO'^Mg 
The corresponding value of stellar mass completeness for i*^ < 24 is given in parenthesis for each bin. 



from redshift z = 0.3 to 1.2. 



(M* > 3 X IO'^Mq, t < 24). By adopting this selection we 
focus on the high-mass end of the galaxy mass function (II- 
bert et al. 2009b). The final samples consist of 28455 opti- 
cally selected galaxies and 8234 LIRGs in the given redshift 
range. 

The redshift distribution of the two samples is shown in Fig- 
ure |3] where the X-ray detected AGN are also highlighted. 
AGN with L(2-10 keV)> 10"*^ erg/s among the 24 ^m detec- 
tions down to 80 /iJy and /"^ < 24 represent a small fraction 
of the population: about 3% in the redshift range z = 0.3 ^ 1 .2 
and 4% at z = 0.8 -^ 1.2. For these the SFR might be over- 
estimated due to the AGN contribution to the infrared SED. 
We will show that the trends are not modified if they are ex- 
cluded in the LIRG sample. We do not attempt to account 
for the contamination by obscured, Compton-thick AGN at 
z ^ 1 . However, we notice that they can hardly outnumber 
X-ray detected AGN (Fiore et al. 2008, Martinez-Sansigre 
et al. 2007). The fraction of highly obscured, not X-ray 
detected AGN at z = 0.7 ^1.2 and with intrinsic luminosity 
L(2- lOkeV) = 10435-44 gjg/g jg 3 7 ^ io5Mpc-\ correspond- 
ing to 67% of that of X-ray selected AGN (Fiore et al. 2009). 

4. PROJECTED DENSITIES 

We adopted the Nth nearest neighbor density estimate 
(Dressier 1980, Capak et al. 2007, Guzzo et al. 2007, Post- 
man et al. 2005, Smith et al. 2005). Projected density mea- 
surements based on photometric redshifts are able to detect 
the structures and reconstruct the overdensities (Scoville et 
al. 2007b). They have been successfully used to detect the 
large scale structure and nicely match the structures detected 
by X-ray and spectroscopic surveys. Although limited by 
the photo-z precision (which implies a smoothing effect), the 
photometric redshifts trace the dependence of galaxy prop- 
erties on environment in the LSS and can use large unbi- 
ased galaxy samples, a factor of ^ 10 larger samples than the 
largest spectroscopic surveys thus limiting the statistical scat- 
ter (Scoville et al. 2007b, Scoville et al. 2010). In particular, 
Scoville et al. (2007b) showed that the correlations between 
galaxy properties, such as SFR, stellar mass, and environment 
are not smoothed out when using accurate photometric red- 
shifts. Conversely, spectroscopic surveys allow for precise 
distance measurements but are limited by the sampling rate 
which is usually low. 

The projected density around each galaxy is derived from 
the distance to the Nth neighbor, r, which defines a circular 
area whose surface density is 'Sn = -t- • This density is then 
divided by the median of projected densities of each redshift 
slice across the whole COSMOS field, in order to obtain a rel- 
ative density Sat/S in each redshift interval. We compute pro- 
jected densities for each galaxy with magnitude i'^ < 24 and 
stellar mass larger than lO'^M©, outside masked areas con- 
taminated by saturated stars (Ilbert et al. 2009a), using 50390 
galaxies in total. We use a cylinder with a depth of ±3a, 
centered on each galaxy. The main difference with previous 



studies on the COSMOS field is that the accuracy of the pho- 
tometric redshifts has now improved by roughly a factor of 3. 
The density field of the COSMOS data is measured using the 
same selection described above. At these stellar masses the 
completeness is 98, 88, 66 and 43% from z =0.3 to 1 .2 (Figure 
121). As discussed also in Scoville et al. (2007b), a higher stel- 
lar mass cut-off would be more conservative and yield larger 
completeness, but at the same time would limit the sample 
to only the most massive galaxies and under-represent lower 
mass galaxies, therefore limiting the dynamic range. There- 
fore when computing densities we allow for a lower complete- 
ness. Conversely, when comparing LIRGs to optical galaxies 
we adopt a cut at M* > 3 x IO'^Mq in order to have a higher 
completeness level (see Sect. 3). 

We used the Millennium Simulation mock catalogs 
(Kitzbichler & White 2007) to compare the derived surface 
densities with the three dimensional (3D) volume densities 
for galaxies selected in the same way as for the COSMOS 
data. We run the simulation here described using the same se- 
lection used for the observed data. First, we compute the 3D 
volume densities by using the positional information avail- 
able from the Millennium Simulation, and measure the den- 
sity field around each galaxy in a sphere of radius equal to the 
distance of the 10th nearest neighbor In the simplest assump- 
tion, i.e. spherical symmetry, the 3D densities are related to 
the projected densities by p = 10 x 3/47rr^(), where rio is the 
distance to the 10th nearest neighbor No redshift errors were 
included in the 3D density computation at this stage. This 
density is related to the projected density by the following 
relation: Sio = ap^^^ where a = 1.76. Figure|4]shows the cor- 
relation between the input and recovered densities. The hori- 
zontal axis shows the projected 3D volume density in the Mil- 
lennium simulation computed as described above. The verti- 
cal axis shows the Eio measured in the Millennium simula- 
tion in the same way as for the COSMOS data as described at 
the beginning of this section. This test indicates that projected 
densities on average overestimate the real densities by a factor 
of 2.2 but there is proportionality between the 2D and the 3D 
density estimates for the simulated data (Figure 4, black con- 
tours). In order to estimate how the photo-z errors affect the 
2D density computation we associate to each redshift in the 
mock catalog a random error drawn from a gaussian distribu- 
tion with (T./(l +z) = 0.012, equal to the error in the measured 
photometric redshift of COSMOS, and we measure the pro- 
jected densities. A gaussian error distribution is an approxi- 
mation of the real error distribution. However, it is reasonable, 
considering that catastrophic errors in the photo-z are not an 
issue both for the optically selected and the 24 fim selected 
samples down to /"^ = 24, since they are of the order 1 % and 
will therefore only affect densities by that amount (Ilbert et al. 
2009a). The effect of the photometric redshift error is shown 
in Figure |4] with red contours. The errors in the photo-z in- 
crease the scatter in the correlation between 3D and projected 
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Fig. 4. — Isodensity contours of measured Sio versus projected 3D galaxy 
density, from the Millennium Simulation (black contours). Red contours: 
same relation where the en'or on photo-z is taken into account. 

densities. However projected densities are still proportional 
to 3D densities and environment can be characterized at least 
in the high density regime. The scatter of the correlation is 
large at small densities but it decreases with increasing den- 
sity. This can affect and smooth out possible existing trends 
by scattering galaxies between contiguous density bins. At 
Eio ^ 10 the scatter is about 100% of the density value, mean- 
ing that a galaxy residing in a region of density 10/Mpc^ will 
be assigned to a density range -^ 20/Mpc^ with a probabil- 
ity of 68.3%. At Eioo >100 the scatter is about 25% of the 
density value. Therefore, for example, a galaxy residing in a 
region of 200/Mpc^, has 68.3% probability of being assigned 
to a density in the range 150 H- 250/Mpc^. Around 50/Mpc^ 
the scatter is ~ 50%. This must be taken into account when 
studying trends as a function of the environment, and the den- 
sity bins must be chosen accordingly. In particular, for the 
four density bins used in this paper (see section 5), the la er- 
ror on the density estimate is 150%, 68%, 42% and 32% of 
the projected density value, for each bin from the low density 
to the high density one, respectively. 

5. RESULTS 
5.1. SFR and SSFR versus environment and redshift 

Figure|5] upper panels, shows the SFR of the LIRGs versus 
the projected density in four redshift bins from z =0.3 to 1.2. 
The increase in the median SFR at z = 1 -^ 1.2 is due to the 
luminosity cut adopted here (L/s> 3 x lO^L©). The increase 
of the specific star-formation rate (SSFR) with redshift is ex- 
pected given the evolution of the stellar mass-SSFR relation 
with redshift (Elbaz et al. 2007). 

The SFR-density relation shows a large scatter and the me- 
dian SFR is constant at all densities. The median SFR per 
LIRGs does not vary with the environment. In our analysis 
we do not find any excess of LIRGs in the high density envi- 
ronments up to z =1. We notice that possible existing trends 
are likely to be erased in the low density regime (first two bins 
in projected densities), but the flat trend observed here is sig- 
nificant in the two highest density bins. This said, we do not 
observe the reversal of the SFR-density relation (i.e. average 



SFR increasing with density, Elbaz et al. 2007, Cooper et al. 
2008) up to z = 1 . At z > 1 the median SFR is compatible 
with a flat dependence on density. This difference may be due 
to the highest accuracy in the density measurements achieved 
by using spectroscopic redshifts. However, the samples of 
Elbaz et al. (2007) and Cooper et al. (2008) are not mass- 
limited and spanned a broader range both in stellar mass and 
SFR, probing down to significantly lower SFR compared to 
our LIRG sample. Therefore, such reversal by might be due to 
the contribution of galaxies with lower stellar mass and SFR. 
Indeed, in the local Universe, the dependence of SFR on en- 
vironment is weaker for massive galaxies (Kauffmann et al. 
2004). For massive LIRGs such reversal, if any, might oc- 
cur at z > 1. Figure |5] lower panels, shows the SSFR versus 
density for each of the previous redshift bins. Here again the 
SSFR is consistent with a constant trend at all densities. 

5.2. LIRG fraction versus environment and redshift 

Figure |6] shows the logarithmic fraction of LIRGs among 
the whole population of galaxies with M* > 3 x 1O'°M0, as 
a function of the projected density from redshift z = 0.3 up to 
z = 1.2, and in two stellar mass bins (M* = 3 x 10'°- 10"Mq 
and M* > 1O"M0). For galaxies at z > 1 the density bins are 
larger than at lower redshift in order to contain a sample large 
enough. Filled circles represent the LIRG fraction among the 
whole opticaUy selected population with M* > 3 x IO'^Mq. 
Open circles represent the same fraction but where the con- 
tribution of X-ray detected AGN has been subtracted from 
the LIRG sample. The error-bars on the measurements were 
computed propagating the poisson errors on the galaxy counts 
and using upper limits in the bins where we have only a few 
counts (Gehrels 1986). The error-bars are large in the highest 
density bins, since these are in general less populated (hosting 
a few tens optical galaxies, and a few LIRGs). 

We find that the LIRG fraction correlates with both redshift 
and environment. /) LIRGs show an overall increase of their 
fraction with increasing redshift. Such rise is predicted by 
the evolution of the LIRG luminosity function(Le Floc'h et 
al. 2005, Caputi et al. 2007, Magnelli et al. 2009) and of 
the mass function by spectral type (Ilbert et al. 2009b). //) 
We find that the LIRG fraction decreases in the high density 
environments. At densities larger than about 80 -^ 100 Mpc"^ 
the LIRG fraction drops compared to the low density envi- 
ronment. The decrease of the LIRG fraction with increas- 
ing density is observed up to z ~ 1.2, but the dependence 
flattens with increasing redshift. The difference between the 
3rd and the 4th density bin decreases from the low redshift 
sample (2.0±0.2 factor at z = 0.3 -f 0.6) to the high redshift 
ones (1.6±0.I6 and I.3±0.I factors at z = 0.7 and 0.9, re- 
spectively), for steflar masses M* = 3 x 10'°- lO^M©. The 
first two density bins should not be considered since in this 
range the uncertainties in the density measurement probably 
erase possible existing trends. 

Figure|6] right panel, the LIRG fraction is shown for galax- 
ies with stellar masses larger than M* = 1O"M0. Again, the 
LIRG fraction decreases towards high densities, and the de- 
crease appears to flatten with increasing redshift. The ra- 
tio between the third and fourth density bin is 2.15±0.2 at 
z = 0.3-^0.6 and it decreases to 1.6 ±0.15 and 1.3 ±0.1 at 
z = 0.7 and 0.9, respectively. At z > 1 the LIRG fraction is 
constant at all densities within the errors for both stellar mass 
bins. 

We note that the dependence is observed also for the LIRG 
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Fig. 5. — Upper panels: SFR versus density for the LIRG sample (black symbols) and AGN (magenta symbols) in redshift bins. Red symbols represent the 
median values and standard eiTors in density bins, omitting AGN. Lower panels: SSFR versus density. 



sample cleaned from contribution by the AGN detected by 
XMM-Newton survey. Indeed, the two samples, with and 
without AGN, follow the same trend within the errors (Figure 
|6|. The environmental properties of AGN in the large-scale 
structure and the comparison with star-forming galaxies are 
beyond the aim of this paper, and they will be addressed in a 
separate work. 

The decreasing fraction observed here is not due to projec- 
tion effects since: a) the size of the two highest density bins 
is larger than the average scatter of the relation shown in Fig- 
ure 4; b) the errors on projected densities tend to smooth out 
trends by scattering galaxies between contiguous bins. We 
explore the possibility that the decreasing fraction of LIRG in 
the high density environment might be produced by mass seg- 
regation, i.e. the most massive galaxies are formed in the most 
massive halos at earlier times, and therefore, preferentially re- 
side in high density environments. As mentioned before, we 
studied the environmental dependence in two stellar mass bins 
in order to minimize such effect. In Figure [T] we plot the nor- 
malized distributions of stellar masses of LIRGs and of the 
optical sample in redshift bins and in two regimes, high and 
low projected density (> and < 90/Mpc^ respectively). The 
red histograms show the LIRG fraction, with poisson errors 
taken into account. Such fraction is plotted in 0.5 dex wide 
bins of stellar mass in order to minimize the noise. The mass 
distribution of LIRGs is not skewed towards lower masses 
compared to optically-selected galaxies. The Kolmogorov- 
Smirnov (KS) test shows that the distributions of LIRGs and 
optical galaxies derive from the same parent population with 
a probability larger than 98% for all the bins except the last 
three (z = 1 -^ 1.2 and z = 0.8 4- 1, high density), for which 
the KS test gives no significant result. However, the LIRG 
fraction does not show any significant trend as a function of 



stellar mass in any environment in the stellar mass range stud- 
ied here. Therefore the decreasing fraction of LIRGs in high 
densities (Figure|6l) is not a mass related effect but it is due to 
the environment. 

Figure [8] shows the fraction of LIRGs as a function of 
redshift, normalized to their values in our first redshift bin 
(z = 0.5). Empty triangles represent the whole sample, re- 
gardless of the environment. The LIRG fraction increases 
with redshift, reflecting the evolution of the whole infrared 
selected population. This increase is similar to the global in- 
crease observed in Figure |6l We compare it with the expec- 
tations from the known evolution of the galaxy population: 
we use the mass functions of Ilbert et al. (2009b) to compute 
the expected number of galaxies with M* > 3 x lO'^M© in 
our field as a function of redshift, and the infrared luminosity 
function of Magnelli et al. (2009) for the number of LIRGs 
in the field. Given the relation between SFR and stellar mass 
at z ^ 1 calibrated by Elbaz et al. (2007), we expect the aver- 
age SFR of a M* = 3 X lO^'M© galaxy will be 20± lOMg/yr. 
With our LIRG selection, equivalent to SFR > 18M0/yr we 
expect that most LIRGs lie above our mass threshold, but not 
all of them, given the scatter of the relation. Therefore we 
can compute from the luminosity and mass functions (LF and 
MF) an upper limit to our observed relation. The solid curve 
in Figure [8] shows the expected upper limit variation for the 
LIRG fraction. We find a similar behavior, with the LIRG 
fraction increasing with redshift, but at a lower rate than what 
the LF and MF would predict: between z = 0.5 and z = 1, 
we observe a 3 fold increase in the LIRG fraction, while the 
computed increase is 4 fold. This discrepancy is expected 
since the solid line is an upper limit. We plot with open cir- 
cles the evolution of the LIRG fraction in low density envi- 
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Fig. 6. — Fraction of LIRGs to the optically selected galaxy population with M, > 3 X IO'^Mq as a function of galaxy density and redshift in two stellar 
mass bins: M, = 3 X 1O"'-^1O"M0 (left panel), M, > 1O"M0 (right panel). Filled circles: X-ray detected AGN included, open circles : X-ray detected AGN 
excluded. 
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Fig. 7. — Normalized distributions of the stellar mass for the optically 
selected sample (thick solid histogram) and for LIRGs (shaded histogram), 
in redshift bins. Left panels show the distributions for projected density 
< 90/Mpc^ , right panels those for density > 90/Mpc^ . Red histograms show 
the LIRG fraction to optically selected galaxies, and poisson errors. 



ronment (5 < Sio < 10 Mpc ), and with solid circles that for 
high density environment (Sio > 90 Mpc"^). We find that the 
LIRG fraction evolves faster (a factor of 2) in high density en- 
vironments than in low density ones, at 2a confidence level. 
However, given the large errors, the high density environment 
fractions are still consistent with the low density environment 
fractions within 3 a. 

6. CONCLUSIONS 

We studied the properties of a sample of LIRGs detected 
at 24 /im in the COSMOS survey and with stellar masses 




Fig. 8. — Evolution of the normalized LIRG fraction with redshift in 3 
redshift bins up to z =1. Open triangles: all LIRGs, filled circles: LIRGs in 
high density environment (Sio > 90), open circles: LIRGs in low density 
environment (Eio = 5 ~ 10, shghtly shifted for clarity), solid line: expected 
evolution from Magnelli et al. (2009) LF and Ilbert et al. (2009b) ME The 
fractions are normalized to the value at z = 0.3 -=- 0.6. 



M* > 3 X 1O'"M0 and z = 0.3 ^ 1 .2, as a function of the en- 
vironment, by using accurate photometric redshifts. In or- 
der to disentangle environment from any other possible ef- 
fect (mainly mass segregation) we studied the properties of 
our infrared galaxy sample in two bins of stellar mass. The 
environment was characterized by using the distance of the 
10th nearest neighbor, and the derived projected galaxy den- 
sity around each galaxy in redshift slices. Our results can be 
summarized as follows: 

• The median SFR and SSFR per LIRG increase with 
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redshift, in agreement with the expected evolution of 
the LIRG luminosity and mass functions (Le Floc'h et 
al. 2005, Caputi et al. 2007). The median SFR and me- 
dian SSFR are constant at all densities within the errors, 
in agreement with what found by Hwang et al. (2010) 
in local LIRGs. We do not observe the reversal of the 
SFR density relation up to z ^ 1 for our mass-selected 
LIRG sample. Such reversal, if any, might occur at 
z > 1 for high-mass LIRGs. This is not necessarily true 
for galaxies of smaller stellar mass. 



• We find that the fraction of massive LIRGs among op- 
tically selected galaxies with the same stellar mass, 
varies as a function of both redshift and environment. 
The global LIRG fraction increases with increasing red- 
shift, as expected from the cosmic evolution of the SFR 
density and of the infrared LF The LIRG fraction de- 
creases with increasing environment density. Such de- 
crease flattens going from z = 0.4 to z ~ 1 and the frac- 
tion is constant in all environments for z > L These 
results imply that the evolution of the LIRG fraction 
among optically selected galaxies from z = 1 to z = 0.4 
is roughly twice as fast in dense than in low density en- 



vironments. 

• Finally, we find that a large fraction of LIRGs is present 
in the most dense environments by z ^ L These might 
be the progenitors of the massive galaxies found today 
in the center of clusters, which therefore were not likely 
formed through dry mergers. 

Using photometric redshift, with their inherent limitations, 
we have shown that LSS and environment has a role on star 
formation. Future works based on spectroscopic surveys will 
refine these results in the densest structures, and open the 
realm of low density where our estimator is limited. 

This work is based on observations made with the Spitzer 
Space Observatory, which is operated by the Jet Propulsion 
Laboratory (JPL) , California Institute of Technology under 
NASA contract 1407. Support for this work was provided 
by ANR "D-SIGALE" ANR-06-BLAN-0170 grant. Sup- 
port for this work was provided in part by NASA through 
contracts 1282612 and 1298231 issued by the JPL, and grants 
1289085, 1310136. We thank the anonymous referee for care- 
fully reading and providing useful suggestions that improved 
the paper 
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